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1. Introduction 

With the Wilson fermion formulation chiral symmetry is explicitly broken at finite lattice spac- 
ing and the consequences of this explicit breaking have to be dealt with. Among those is the fact 
that 0(a) counter terms in the Symanzik expansion are no longer excluded and that there is no 
conserved Noether current associated with a continuous chiral symmetry of the lattice action. 

The first problem is addressed by the Symanzik improvement programme, i.e. the addition 
of irrelevant operators to both the action and the composite fields. By tuning the coefficients of 
only a few terms it is possible to remove the 0(a) scaling violations from the theory. Here we 
compute the improvement coefficient for the isovector axial current ca(^o), which - together with 
the improvement of the action [1] - will ensure the absence of 0(a) lattice artifacts in all matrix 
elements of the PC AC relation, which involve insertions at finite separation [2]. 

The normalization factor Za(#q) of the improved axial current is obtained by enforcing a 
continuum-like transformation behavior under infinitesimal chiral rotations. Since the isovector 
chiral symmetry is recovered in the continuum and the normalization condition is based on a local 
identity, Z\ is finite and scale-independent. 

For unexplained notation and additional details about the improvement and normalization con- 
ditions we refer to [3] and [4], respectively. 

2. Axial current improvement 

When calculating a bare quark mass on the lattice from matrix elements of the PCAC relation 

dpAfa) =2mP a (x) , (2.1) 

any dependence on the kinematic parameters or external states, i.e. the dependence on the precise 
choice of matrix element, is a cutoff effect. A non-perturbative improvement condition can thus be 
obtained by inserting into eq. (EtJ) the improved axial current [2] 



(Ai)£(jc) =A a tl {x) + ac K \(d tl + d;)P a {x) , (2.2) 



and tuning ca such that the masses obtained from two different matrix elements agree. In ( |2.2| ) d„ 
(<?*) denotes the forward (backward) lattice derivative. 

When evaluating improvement coefficients non-perturbatively one has to keep in mind that 
due to cutoff effects in the correlation function used to formulate the improvement condition, the 
coefficients themselves are uncertain to 0(a). 1 While this forbids a unique definition of the im- 
proved theory, the 0(a) ambiguities can be made to disappear smoothly if the improvement con- 
dition is evaluated with all physical scales kept fixed, while only the lattice spacing is varied [5]. 
In the evaluation of our improvement condition we keep the bare quark mass (using the 1-loop 
value for ca from [6]) constant, thus ignoring small changes in renormalization factors, and fix 
the relative lattice spacing in the range of bare gauge couplings we consider through asymptotic 
scaling [3]. 

It is also important to make sure that the correlation functions in the improvement condition 
are not dominated by states with energy close to the cutoff. If this were the case, the improvement 



'For the same reason is uncertain to 0(a 2 ) after improvement. 
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condition might cancel exceptionally large scaling violations and a c A obtained in this way could 
introduce significant 0{a 2 ) effects. 



We construct matrix elements of (2.1) between pseudo-scalar states and the vacuum in the 



Schrodinger functional [7, 8]. More precisely, in this work we use 

f A (x ;G>) = -|^I> (*)^(co)) , (2.3) 
and /p(x ;t») = -^Kn^W) , (2-4) 

with the pseudo-scalar operator 

^( W )=a 6 £C(x)7 5 T c '^(x-y)C(y). (2.5) 

It lives at the xq = boundary of the SF cylinder and depends on a spatial trial "wave function" CO. 
A quark mass from ( |2~l| ) using ( [2J1 ) is then given by m = r + ac A s + 0(a 2 ) with 

\{d {) + d*)f A {x ;co) d dSMxo;(o) 

K*o;<») = - TTT-, — —s and s{x Q ;co ) = — —, — —— ■ (2.6) 

2/ P (^o;ft>) 2/ P (^o;ft>) 

In our chosen setup we now have xq, the insertion time, and CO, the spatial trial wave function, as 

parameters for probing the PCAC relation. Enforcing the independence of the quark mass from 

these results in . , , . , . 

Ar 1 r(%o;fiVi) -r(y \CD n fl) 

-ca = — 7- = t \ - t r , (2.7) 

aAs a s(aco;fi) ff ,i)-jO'o;fi%t ) o) 

and therefore the sensitivity to c A is given by aAs. 

The combined requirement of large sensitivity to c\ and explicit control over excited-state 
contributions is met with the method proposed in [9] where it is tested in the quenched approxi- 



mation. There CO^o/i from (g/7J) are chosen such that &{(Q K $/\) excite mostly the ground and first 
excited state in the pseudo-scalar channel, respectively. Thus the energy of the states can be moni- 
tored directly using the effective mass of / A / P and the sensitivity aAs = as(xo, O^i) — as(xo, &Vo) x 
m* 2 — nv\ is also expected to be large. 

The simulations listed in Table |l| are at constant physical volume according to asymptotic 
scaling as described in [3] and the bare (unimproved) PCAC mass is kept constant to ~ 10%. 

We simulate a set of three spatial wave functions and 
approximate the combinations ^(o w /i) that project to 
the ground and first excited states through the eigenvectors 
of the correlation matrix [3,9] 



L/a 


P 


m L 


-CA 


12 


5.20 


0.18(1) 


0.0638(23) 


16 


5.42 


0.200(5) 


0.0420(21) 


24 


5.70 


0.182(5) 


0.0243(36) 



Table 1: Simulation parameters for c A . Mco',CO) = —L(^(co') a (co)) , (2.8) 

where C is a pseudo-scalar operator at the xo = T boundary. In Fig. |l| two distinct signals are 
clearly visible, which indicates that the approximate projection method works well at our parame- 
ters. The energy of the first excited state is not far away from aT l , suggesting that in even smaller 
volumes the residual 0(a 2 ) effects would grow rapidly. 



Our definition of c A is completed by fixing x$ = yo = T/2 in (2.7) and specifying L = T and 



6=0. The results from the simulations summarized in Table [T] are shown in Fig. |as a function 
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CO 




1.00 




Figure 1: Effective pseudo-scalar masses of fp Figure 2: Simulation results for c^. The solid 
from the j5 =5.42 run. 



line is given by the interpolating formula (2.9) 



of gg. The solid line is a smooth interpolation of the simulation data, constrained in addition by 
1-loop perturbation theory: 

1-0.4485^ 



ca(So) 



-0.00756sjj x 



(2.9) 



1-0.8098^ • 

It is our final result, valid in the range 0.98 <gl < 1.16 within the errors of the data points (at most 
0.004). By performing additional simulations we have verified that the volumes in our runs are 
scaled sufficiently precisely such that systematic errors due to deviations from the constant physics 
condition can be neglected. The same also applies to variations in the quark mass. 



3. Axial current renormalization 

In a massless renormalization scheme preservation of 0(a)-improvement implies that the 
renormalized improved current is of the form [2, 10] 

(A R )£=Z A (l+ft A m q )(A T )£. (3.1) 

The normalization condition we use [4, 11] is based on [12], the ALPHA collaboration's quenched 
determination of Z A . Since a massless scheme requires the normalization condition to be set up at 
vanishing quark mass, in [12] the mass term of the axial Ward identity was dropped in the derivation 
of the normalization condition. In practice this condition shows a very pronounced quark mass 
dependence and thus a potential extrapolation is rather steep and the location of the critical point 
must be known with high precision. 

Performing a chiral transformation in the continuum and keeping track of the mass term results 
in the integrated Ward identity [4] 



jd 3 yd 3 xe abc ^A a (y +t,x)A b (y)^ t 



-2m ld 3 yd 3 xr + L Q£ abc (p a (x)A b (y)0 ext ) = i /d 3 y ( \' ( '/( v)^ cxl 



(3.2) 



and the choice 2 



-E cde €>" d & e /6L 6 allows us to replace the r.h.s. of (3.2) with f\ by using 



isospin symmetry. A normalization condition is then obtained by requiring that (3.2) holds for the 



renormalized improved current (3.1 ). 



"No use of spatial wave functions is made in the axial current renormalization, i.e. here i 



ft) = const 
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Figure 3: Comparison of the chiral extrapola- 
tion at j3 =5.2 using the new and old normaliza- 
tion conditions forZA- 



Figure 4: Result for Za from numerical simu- 
lations and 1-loop perturbation theory (dashed 



line). The Pade fit (solid lines) is given by (3.3). 



Fig. [| shows the results of an evaluation of this condition at a bare gauge coupling of j8 = 5.2. 
To show the effect of the inclusion of the mass term, also the results with the method from [12] are 



shown, i.e. when dropping the second term on the l.h.s. of ( |3.2| ). While for the new normalization 
condition the slope in am is consistent with zero 3 , the estimate of Za from the old condition changes 
by 30% in the (small) mass range shown. We anyway see that for am < 0.02 all mass effects show 
a linear behavior. For j3 = 5.5 the extrapolation is similar to the one shown and at all other gauge 
couplings we can in fact interpolate from two simulations very close to the critical point. 

In addition to three (again asymptoti- 
cally) matched lattice sizes at L/a = 8, 12 
and 16, we simulated at three larger values 
of j6 and fixed L/a = 8, which corresponds 
to very small volumes. This was done in 
order to verify that our non-perturbative 
estimate smoothly connects to the pertur- 
bative predictions [13, 14]. The first er- 
ror we quote for Za is statistical and the 
second represents our estimate of the sys- 
tematic error, which originates from devi- 
ations from the constant physics condition. There is no estimate of the systematic error for the 
/3 =5.29 run, which was done only to verify the rapid change of Za in this region of bare gauge 
coupling. It is thus also excluded from a fit, which results in the interpolating formula (again 
incorporating the 1-loop asymptotic constraint) 

„ 2 1-0.918^ + 0.062^ + 0.020^ 



P 


L/a 


T/a 




Z A 


5.200 


8 


18 


0.135856(18) 


0.7141(123) 


5.500 


12 


27 


0.136733(8) 


0.7882(35)(39) 


5.715 


16 


36 


0.136688(11) 


0.8037(38)(7) 


5.290 


8 


18 


0.136310(22) 


0.7532(79) 


7.200 


8 


18 


0.134220(21) 


0.8702(16)(7) 


8.400 


8 


18 


0.132584(7) 


0.899 1(25X7) 


9.600 


8 


18 


0.131405(3) 


0.9132(H)(7) 



Table 2: Results for the chiral extrapolations of Za and 
estimates for the critical hopping parameter K c . 



1 



0.8015 



(3.3) 



4. Summary and Outlook 



For the O(a) -improved action with non-perturbative csw [1]> we have determined the improve- 
ment coefficient ca for j8 > 5.2, which roughly corresponds to a < 0.1 fm. The improvement condi- 
tion was implemented at constant physics, which is necessary in a situation when O(a) ambiguities 
in the improvement coefficients are not negligible. 



Any remnant slope is due to the neglect of as well as contact terms in the second term in the l.h.s. of (3.2) 
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Through the calculation of Za(#o) we have shown that in a lattice theory with two flavors of 
Wilson fermions normalization conditions can be imposed at the non-perturbative level such that 
isovector chiral symmetries are realized in the continuum limit. Since we are working with an 
improved theory, chiral Ward-Takahashi identities are then satisfied up to 0(a 2 ) at finite lattice 
spacing. 

The improvement and normalization conditions were implemented in terms of correlation 
functions in the Schrodinger functional framework and evaluated on a line of constant physics 
in order to achieve a smooth disappearance of the 0(a) and 0(a 2 ) uncertainties. Clearly, the meth- 
ods employed in this paper may also be useful to compute ca and Za in the three flavor case, where 
csw i s known non-perturbatively with plaquette and Iwasaki gauge actions [15-17]. 

The determinations of ca and Za were carried out within the ALPHA collaboration's pro- 
gramme to calculate quark masses in a fully non-perturbative framework. The results of this pro- 
gramme and its application to the strange quark mass are presented in [18, 19]. 
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